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Abstract. We study graphene nanoribbons (GNRs) with armchair edges in the 
presence of Rashba spin-orbit interaction (RSOI). We impose the boundary con- 
ditions on the tight binding Hamiltonians for bulk graphene with RSOI by means 
of a sine transform and use the results to derive different approximations and to 
investigate their ranges of validity. Finally, we derive an approximation for the 
lowest two energy bands, which is valid for experimentally available sizes of RSOI. 
In addition we study the spin polarization in GNRs resulting from the RSOI. 
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1. Introduction 

Spintronics is a multidisciplinary field whose central theme is the active manipulation 
of the spin degree-of- freedom in solid state systems [1]. The integration of spintronics 
concepts into standard electronic devices is an important technological challenge for 
the semiconductor industry. The goal is the realization of devices, which for equal 
sizes, have better performance and lower power consumption with respect to state- 
of-the-art electronic devices [2]. In this respect, graphene has great potential for 
spintronics applications [3]. On the one hand it is characterized by very weak spin- 
orbit interactions (SOIs) [4] making pristine graphene the spintronics material with 
the longest spin-coherence time [5]. On the other hand, there are several promising 
proposals to manipulate the spin states via local material engineering of the graphene 
membrane [6, 7]. 

There are two possible types of SOIs in graphene: intrinsic and Rashba (R). Both 
of them can be understood in terms of the symmetry properties of the honeycomb 
lattice [8] and by tight binding arguments [4]. The intrinsic SOI opens a gap in 
the energy spectrum and transforms graphene into a two-dimensional topological 
insulator [8]. However, intrinsic SOI is extremely weak in prestine graphene. Armchair 
GNRs with intrinsic SOI have been studied in reference [9]. 

The RSOI is mainly connected to the overlap of the it and the a orbitals of 
the carbon atoms in the sp 2 -hybridization. It can be manipulated by an electric 
field perpendicular to the graphene plain or by the local curvature of the graphene 
sheet [4]. Alternatively, it was proposed that RSOI can be enhanced by modifying 
the bonding within the graphene sheet by covering it with, e.g., hydrogen [10]. It 
has been predicted that the first two methods lead to a RSOI strength of A ~ lmeV 
while the third method could lead to a RSOI strength of approximately one order of 
magnitude larger. In an experimental study [7] it was shown that opportune substrate 
engineering could even lead to A ~ 200 meV. Thus there is quiet some reason to hope 
that manipulation of RSOI in graphene will soon become good enough to realize RSOI 
based spintronics devices. In this respect several efforts have already been made in 
order to investigate the spintronics properties of bulk graphene, e.g. in the case of 
inhomogeneous SOI structures [11, 12, 13, 14, 15]. 

However, particularly interesting for graphene is the stripe geometry, the so-called 
graphene nanoribbon (GNR). While large graphene sheets are still polycrystalline 
and far from perfect, GNRs with an armchair edge (cf. figure 1) have already been 
synthesized with atomic precision and well defined edges [16]. In this case, the dangling 
a-bonds of the GNRs were passivated with hydrogen, so that also at the edges the 
sp 2 -hybridization is preserved. 

So far GNRs with RSOI and zig-zag edges have been intensely studied. For 
example, Zarea and Sandler [17] derived the analytic energy spectrum of a zig-zag 
GNR with RSOI produced by an electric field perpendicular to the GNR surface. 
Gosalbez-Martmez et al. [18] studied RSOI produced by effects of curvature in zig- 
zag GNRs. Also GNRs with hardwall boundary conditions and RSOI received some 
attention [19]. On the contrary, armchair GNRs with RSOI received little attention 
so far — despite of their experimental availability [16]. 

In this article, we fill the gap and study armchair GNRs with RSOI using an exact 
and different approximative methods. The article is organized as follows. Section 2 
presents the system studied and reviews the ingredients needed to derive the RSOI 
Hamiltonian for an armchair GNR in section 3.1. We linearize the nearest-neighbor 
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Figure 1. An armchair GNR of width W = 4a is shown, where a is the width of 
the hexagon. The green and red dots mark the A and B atoms of the honeycomb 
lattice, respectively. Here d\ , d,2 and ^3 are the displacement vectors connecting 
nearest neighbor sites. The yellow square indicates the unit cell of the ribbon, 
with 4x4 + 2 atoms for the width shown. The distance between the first missing 
rows of carbon atoms is denoted by W = W + a. The inset depicts the first 
Brillouin zone of bulk graphene with the two inequivalcnt K, K' points. 



hopping Hamiltonian and the RSOI Hamiltonian in order to investigate armchair 
GNRs in the Dirac approximation in section 3.2 and obtain an analytic equation for 
the lowest energy bands in section 3.3. Finally, in section 3.4 we investigate the spin 
polarization arising due to the RSOI in the ribbon. 

2. Basic Theory 

2.1. Graphene armchair nanoribbons 

We study armchair GNRs with the edge aligned along the y-dircction as depicted in 
figure I. Here the two inequivalent carbon atoms A and B are marked by green and 
red dots, respectively. The three displacement vectors which connect nearest neighbor 
sites are 

Qj , . . a I v3 1 \ 

*-^(m). ^=^{±-,- 2 y a) 

where a = ^/2la cc is the width of the hexagon and a cc = 0.142 nm the carbon-carbon 
distance. For all plots we choose a cc as our unit of length and formally set a cc = 1. 

The first Brillouin zone of bulk graphene, which is also a hexagon rotated by tt/2 
with respect to the graphene lattice, is shown in the inset of figure 1. The wavevectors 
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in the x and y direction are q x and q y , respectively. There are two inequivalent 
Dirac points, called K and K', at which the conduction and valence bands touch 
and, in the absence of RSOI, disperse in a conical shape. For the expansion of the 
spectrum around the Dirac points, we will use the two points on the q x axis, namely 
K = (47r/(3a), 0) = (K x , 0) and K' = -K. 

There are four atoms for every column of the ribbon and two additional atoms in 
order to truncate the ribbon in an armchair configuration. Note that in this article, 
the distance between the first missing atoms W will be used instead of the actual 
width of the ribbon which is W — W — a. 

2.2. Tight binding Hamiltonians 

The Hamiltonian in the tight binding approximation is the starting point of our 
treatment. The kinetic nearest-neighbor hopping Hamiltonian is given by [20] 

n kin = - t J2 X>1 (-Hi + djK(-H<) + H.c] (2) 

i,j a 

with the hopping integral t w 2.7 eV. The operator a a (Ri) annihilates a quasiparticlc 
on the A atom at lattice position Ri and with spin a and bURi + dj) creates a 
quasiparticle at the B atom at position (_R,; +dj) with spin a. The sum over j = 1, 2, 3 
takes into account the three displacement vectors. 

The RSOI Hamiltonian induced by, e.g. an electric field perpendicular to the 
graphene sheet, is given by [17, 21] 

H so = iXj2J2 b t'( R i + d A 

i,J cr.cr' 

and has the form of a spin-dependent nearest-neighbor hopping with hopping integral 
A. Here z is the unit vector in ^-direction, s is the vector of Pauli matrices associated 
with the spin degree of freedom and dj are the normalized displacement vectors. The 
Hamiltonian (3) lifts spin degeneracy, since hopping between the A and B atoms is 
always accompanied by a spin-flip. 

2.3. The Dirac approximation 

For the Dirac approximation of the nearest-neighbor hopping Hamiltonian (2) the 
Fourier transform is taken and expanded close to the Dirac points K and K' . Thus, 
q x = ±K X + k x and k y = q y , where k x and k y are small compared to \K\ [20]. This 
leads to 

^Dirac = V F (T z S a x k x - T Q S (T y k y ), (4) 

where vp = (y/3at)/2 is the Fermi velocity, o~i, Sj and tj (i = x, y, z) are Pauli matrices 
representing the A/B-sublattices, the fermionic spin, and the K / K 1 valley degree of 
freedom, respectively. The index is used for the unit matrix. Note that %Dirac 
contains two blocks 'Hp irac and 'H§ irac related by time- reversal symmetry (TRS): 
^Dtac(^) = ^Dirac( — This facilitates the calculation of the spectrum for an 
armchair GNR, as we will show in 2.4. 

In order to point out some difficulties that arise when RSOI in the Dirac 
approximation is considered together with armchair boundary conditions, we give the 
RSOI Hamiltonian up to zeroth order in the wavevector fe (see [4, 8] and Appendix A): 
3A 

^so = -7r{ T oo-xSy - T z a y s x ). (5) 



a^R^ + R.c. (3) 




Figure 2. Figure (a) shows the energy dispersion of bulk graphene in the presence 
of RSOI in dependence of the wavevectors q x and q y . The cut along the <j x -axis, 
marked in red, is shown in figure (b). The blue, dashed line is calculated taking 
into account Wso U P to first order in k, for the green, dashed line only the zeroth 
order in k is used for Hso- Figure (c)-(e) show the same cut for different RSOI 
strength, A. 



As before, the Hamiltonian obeys TRS. But now it explicitly depends on the spin. 
Therefore it is impossible to obtain a quantization condition for the wavevectors of a 
ribbon with RSOI in the low energy approximation, like it is done in section 2.4 for 
the Dirac Hamiltonian. Therefore, in section 3, we will start from the tight binding 
Hamiltonians and instead of taking a Fourier transform, we will take a sine transform 
to ensure that the boundary conditions are fulfilled. Afterwards the linearization can 
be done without problems. 

Before we proceed, let us focus on RSOI in bulk graphene for later comparison. 
The energy spectrum of bulk graphene for different A is shown in figure 2. The spectra 
shown in red are derived from the tight binding Hamiltonians (2) and (3), for the 
approximations to the spectra shown in light green, the linearized nearest-neighbor 
hopping Hamiltonian (4) and the zeroth order RSOI Hamiltonian (5) are used. One 
can see, that the light green spectra do not catch the triangular warping which is 
observed for A ~ 0.5< and larger. In this case, taking the linear order expansion of the 
RSOI Hamiltonian (blue, dashed lines) works slightly better. For a derivation of the 
linear order, see Appendix A or reference [22]. 
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2-4- Armchair boundary conditions for GNRs without RSOI 

The boundary conditions for armchair GNRs without RSOI in the Dirac 
approximation where first studied by Brey and Fertig [23]. We briefly review their 
results. 

The wave function on the graphene lattice for energies close to the Dirac points can 
be written in terms of slowly varying envelope functions <fi K and 4> K . Without taking 
into account the spin degree of freedom, these functions are two-component spinors 
and eigenstates of one of the two equivalent upper blocks or one of the two lower 
blocks of equation (4) in the (aK,bx) or (clk 1 ; bx 1 ) representation (cf. Appendix A 
and reference [20]). The wave function on sub-lattice \i = A, B can be written as 

^{x,y) =4 K '*tf(x,v)+e- lK '*4>F (x,y). (6) 
In this approximation the variables x and y arc continuous. The two functions cf>^ 
and 4>^ are expansions of the full wave function at the two different Dirac points K 
and K 1 and the fast varying phases e xX are required in order to bring both partial 
wave functions back into the same coordinate system (q x = ±K + k x ). 

For a ribbon that is cut at x = and x = W, all wave functions on the missing 
atoms have to be zero (cf. figure 1). This leads to the boundary conditions 

Mo,v) =<(o ! y) + <'(o ! y) = o, (7) 

MW,V) = e iK * w 4>«(W,y) + e~ iK * w 4>f (W,y) = 0. (8) 
The boundary condition (7) requires that spinors propagating along the x-axis from 
the K valley are superimposed with spinors propagating in the opposite direction from 
the K' valley and vise versa, since the corresponding eigenvectors are equal due to 
TRS. The boundary condition (8) gives the quantization condition on k x , 

fcW = k - ±K + — - — ± (9) 
k x _ kn -±K x +^- w ± (J) 

The integer m labels different transverse modes going from m = 1 to m max , which is 
the total number of modes in a ribbon. The total number of modes, on the other hand, 
is equal to half the number of atoms in a unit cell, indicated by the yellow square in 
figure 1. For the second equality sign we rewrite W = a(3j+7) with a positive integer j 
and 7 = — 1, or 1. The modes n range from — Aj to n max = m max —£j = 4(2 j +7) — 2. 
The shifted modes n are convenient to work with, because n = corresponds to the 
lowest energy bands. The equation for the energy spectrum of the armchair GNR now 



reads E n , t = eyk^ + k\. Only for 7 = the GNR is metallic, since a zero-energy 

solution is possible for k y = 0. For 7 = ±1, the spectrum shows an energy gap equal 
to 2|gi| = 2tt/W. In the following we will focus on the effect of RSOI on metallic 
GNRs, though the method employed is general and also applicable if 7 7^ 0. 

Finally using the q m as wavevector the eigenstates of the kinetic Hamiltonian that 
obey the boundary conditions (7) and (8) are proportional to sine functions, 

* e ,„(x,y) oc e ik y y ^ !kyte sm(q m x) . (10) 

Here, $^ k are the two-component eigenvectors of the Dirac Hamiltonian of the K 
valley, depending on the wavevectors of the exponential ansatz and on the sign of the 
eigenenergy e. The exact form of these vectors is not needed because we only need 
to know that the boundary conditions arc fulfilled by the sine function (10), to give a 
physical interpretation to the general sine transform introduced in the next section. 
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Figure 3. Spectra for two different widths and A = 0.5t are shown. The exact 
spectra obtained with the sine transform are shown in red. The dark green 
spectra are calculated using the continuum approximation, equation (16). The 
two smallest possible widths of metallic GNRs are shown while the spectrum for 
the same A and W = 12a can be found in figure 4 (b). 



3. Armchair GNRs with RSOI 

3.1. Sine transform and continuum limit 

It is important to recall here that the tight binding Hamiltonians are expanded over a 
basis set of counter propagating exponential wave functions to obtain (5) and (4). By 
definition these states do not fulfill the armchair boundary conditions. For this reason, 
in order to study the effect of RSOI for the case of an armchair GNR we need to expand 
the Hamiltonians on a basis set that is naturally fulfilling the boundary condition. 
Therefore we take the sine transform of the a CT (J2j) and b a {Ri + dj) operators in the 
tight binding Hamiltonians (2) and (3). Taking the sine transform corresponds to 
choosing a set of basis functions oc sin(q m x) with the mode q m = nm/W . We use 



/OO I 2 m max 
dq y e lqvm J ^ sm(q m x l )a r7 (q m ,q y ), (11) 

-oo V m max m= i 

/oo I 7\ "W t 
dq y e^ +d ^\ Y, Mdmixi+d^)) 

-oo V m max m —i 

xb a (q m ,q y ), (12) 

where Ri = (xi,yi) and dj = (dj^ x ,dj_ y ). The ribbon is assumed to be infinite 
in y-direction, therefore a continuous wavevector q y can be used. Inserting the 
transformation (12) into equation (3) and using q m = rmr/W yields 

n so = \J dq y Y E {(^' + tf-Mm,q y )a a {q m .,q y ) 
+ + * ai(q m ,q v )b^(q m ,,q y )} (13) 

A™% = ie- iq y (l - e 1 ^ 2 cos(V3g m /2)) S n , m >, (14a) 

/ 1 — (— 1 \ m + m ' \ 

A2% = V3e^/ 2 sin(^3g m /2)£ ro , m , , (146) 



with 
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where 8 mm i is the Kronccker delta and 

1 / 7r(m + ml) \ ^ ^ / 7r(m — m') , 

™' m ™ max + l|_ \ 2 ( m ma X + l)/ \2(m max + l)/J 

We recall that throughout the paper the carbon-carbon-distance a cc is chosen as 
the unit of length. In the limit of a very wide ribbon we can consider continuous 
coordinates x and y and use 

4 to! 

lim B mm i = - -. (16) 

m m „->oo 7r [my — m z 

There arc two terms in the RSOI Hamiltonian: an intra-mode and an inter-mode 
coupling term. Both couple states with different spin, however, the former is coupling 
modes with the same mode index m = m', while the latter couples different modes 
m =/= m! with different parity — i.e. odd modes with even modes. The intra-mode 
coupling is a peculiarity of GNRs associated to the two inequivalcnt carbon atoms in 
the unit cell. No similar effect is present in quantum wires with RSOI [24, 25]. 

Inserting the transformation (11) and (12) in 'Hkin, leads to the well known 
expression: 

/OO m mix 
dq y J2 c~ ia iy/^ (17) 
-°° m=l 



2e iV W 2 cos (^i)j bl{q m ,q y )a a {q m ,q y ) 



where, of course, there is no coupling between different modes. The coupling of the 
two inequivalent valleys K/K' is introduced only by the boundary conditions. 

In the figures 3 and 4 the exact spectra are shown in red and the spectra 
obtained with the continuum approximation (16) in green for different widths and 
RSOI strengths. In figure 3 the two smallest ribbons are shown with A = 0.5t. For 
the smallest conducting width W = 3a the agreement is not perfect, but for W — 6a 
the spectrum shows only small deviations for some of the bands. The magnitude of 
the error introduced by the continuum approximation depends on the value of A as 
illustrated in figure 4. The smaller A, the better the agreement with the continuum 
approximation, because the case without RSOI has no inter-mode coupling and thus 
the continuum approximation corresponds to the exact case. 

We note in passing that it is remarkable how strongly the shape of the spectra 
depends on the widths for strong RSOI coupling. We observe three (figure 4 (b)), two 
(figure 3 (a)) and no Dirac points (figure 3 (b)) for the same A = 0.5 and different 
width. It can even be the case that RSOI opens a gap in the spectrum as in figures 3 
(b) and 4 (c). 

In the next section we will linearize the Hamiltonians (17) and (13) and investigate 
for which parameter ranges such a simplified model is valid. 



3.2. Linearized "Hso an d Hkin for armchair GNR 

The linearized version of both the nearest-neighbor hopping Hamiltonian obtained in 
equation (17) and the RSOI Hamiltonian (13), have the symmetry H^ in / SO (k n ) = 

^kin/so(~^")- Therefore we will only give the linearized versions for the K point 
with k n and k y being the small displacements from the K point, i.e. q m — K x + k n 
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Figure 4. Spectra for different A and W = 12a are shown. The tight 
binding results are depicted by full red lines and compared to the continuum 
approximation using equation (16), depicted by dark green, dashed lines. 



and q y = K y + k y . A straight forward Taylor series of (17) leads to the well-known 
Dirac Hamiltonian with quantized wavevector, 

[%kin]„ = v F s (o- x k n -a y k y ), (18) 

with "H kin = / dq y J2 n [^kin]„ if> n and ip n = {a t (fc„, k y ), 6-f 

(kn ; ky ) ; Oj], (^n ; ky ) , (k n 1 ky 

On the other hand, linearizing the coupling matrices (14a) and (14&) around the K- 
point yields 

0~m,m' ( ~W H ) 




7r(m' — to) 



k?n km' I"' 



Now, the linearized version of the intra-mode RSOI Hamiltonian (13) for modes n = n' 
reads 

„ = y s x [a y + ^ + ^ fc " ) (20 ) 

with "Hg n o ra = / rfli/En'/'n [^so ra ]„ 1 / , '»- Comparing the constant terms to the 
linearized RSOI-Hamiltonian of the infinite graphene plane T-L%q, cf. equation (5), 
we observe that this is the same as adding the blocks for the K and the K' valley 
and dividing by 1/2. This illustrates the coupling of the two valleys induced by the 
armchair boundary conditions. 
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Figure 5. In the upper panel, figure (a)-(c), spectra for different width and 
A = O.li arc shown. Figure (d), in the lower panel, shows the spectrum for 
A = 0.5i and W = 12a. The red, full lines depict the exact result while the blue 
lines show the results after taking the continuum limit and linearizing. The effect 
of the brake down of the continuum approximation is quite pronounced, as shown 
in figure (a), where only the lowest energy bands agree with the tight binding 
result. Also for larger ribbons, as shown in figure (b) and (c), the linearization 
leads to discrepancies. Nevertheless, in figure (c) the bands connected to mode 
and 1 are already in a good agreement with the tight binding method. For 
large A, such as shown in figure (d), the results obtained after linearizing the 
sine transformed Hamiltonians differ substantially from the tight binding result 
especially in comparison with figure 4(b) 



For small but different k n ^ k' n and small k y the elements of the inter-mode RSOI 
Hamiltonian are given by 



["^sc 



(21) 



_ 3 A(1 -(-!)"' + Oyky _ Q- X (k n + k nl ) 

2n\{n' — n) 

where W^q 1 = J dq y J2 n n' Wi [^so° r ]„ „/ Vv- The main and surprising point about 
equation (21) is that there is a contribution coupling even and odd bands that does 
not dependent on W but only scales inversely with the distance of the bands. 
In figure 5 we show the dependence on the width W for a constant A = O.lt. We 
have to use smaller A than in the previous section, since the agreement for A = 0.5i 
is already quite bad, as shown in figure 5 (d). Although we have not made any 
additional restriction on A compared to the preceding section, it is clear from figure 2, 
which shows spectra for bulk graphene, that, as soon as we come in the regime of 
triangular warping, the linear approximation to RSOI also catches the behavior only 
qualitatively. Thus, the linearization is only useful for the low energy and small A 
regime. Since up to now only very small RSOI strengths of about a hundredth of t 
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Figure 6. Spectra obtained with equation (22) are shown in light green, together 
with results of the tight binding method for different width. The red lines 
correspond to W = 12a, the lines colored in magenta correspond to W = 6a 
and the violet lines to W = 3a. 



arc experimentally available, we will go even further in the next section and derive an 
analytic approximation for the lowest energy bands only. 

3.3. Analytic approximation for the lowest bands 

In physical systems the achieved spin-orbit interactions are approximately two orders 
of magnitude smaller as the hopping integral [7]. Therefore we want to propose a 
simplified model for the lowest mode bands (n = 0). Since the coupling between 
the bands is as strong as the coupling between spins, which can be seen in equations 
(21) and (20), and since the band corresponding to n — does couple to the bands 
corresponding to n — 1 this model is only valid far away from band crossings. Adding 
[%km]„ and ["%so ra ] f° r n = and diagonalizing leads to: 

o / k 2 \ 2 
E = ±-^(k y ±\)* + ^ r , (22) 

where we have assumed that t = 1, a cc = 1. Equation (22) shows that the main 
correction due to RSOI is the lifting of the spin degeneracy for the n = mode bands 
and a splitting of the Dirac point into two. Furthermore, there is a term proportional 
to kyX 2 , which opens a gap in the spectrum. But since in the current approximation 
A and k y are both small, this term will only have a neglectable effect. Note also, 
that because we consider the metallic mode (n = 0), equation (22) is independent 
of W. This is demonstrated in figure 6, where equation (22) (light green, dashed 
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lines) as well as the result of the tight binding method for different W (different 
tints of red) is shown. For both RSOI strengths that arc considered, the low energy 
behavior is captured reasonably by equation (22), but the figure also indicates that, 
if A > 0.05i, a better approximation, like the one in the preceding section, should 
be used. Interestingly, for such small A the width dependence for the modes n > 
only shows up in the dependence of k n on W. Therefore, for W = 12a the bands 
corresponding to the modes n = 0, 1, 2.. arc visible in the energy range shown, while 
for W — 6a only half the modes are visible and for W = 3a only one quarter. The 
mode n = 2 for W = 12a gives the same bands as the mode n = 1 for W = 6a. This 
indicates that the width dependence comes from k n only, in contrast to, e.g. figure 3, 
where A = 0.5t and where the width affects the spectra strongly. 

3.4- Spin polarization 

Here, we study the spin polarization, i.e., the expectation value of the x-component of 
the spin operator for different bands [(S x = ^s z ®<Jo)] n , as a function of the momentum 
along the GNR with the tight binding method. 

In figures 7 (a) and (b) we show the results for two different values of RSOI, 
A ~ O.Olf and A ~ O.lt. Similar to the observations for quantum wires in presence 
of RSOI [24, 25], the spin polarization shows a smooth crossover to the opposite sign 
when changing q y from positive to negative values. In the case of the lowest energy 
band, shown in red (n = 0), the crossover is due to a crossing with a band of the 
next higher mode (n — 1), which is lifted by the terms coupling different modes in 
"Hso- The violet and the magenta bands have the same mode n = 1, but still they 
show a smooth crossover from positive to negative spin polarization, which cannot be 
directly attributed to avoided band crossings. Indeed, in the case of quantum wires 
with RSOI [24, 25] it is possible to split the system Hamiltonian in a part with RSOI 
along the wire direction and another one perpendicular to the first. In this respect, 
it is possible to interpret the effects of spin polarization as a function of the mode 
coupling due to the second term. This separation of the Hamiltonian in different 
subparts is not possible in the case of a Dirac-likc Hamiltonian. However, in analogy 
to the case of quantum wires, we observe that the spin polarization curves do not have 
a well defined spin quantization axis. For smaller A the situation is different as then 
the influence of the higher mode bands is negligible, since the modes are energetically 
far apart. This is shown in figure 7 (a). The crossing point of the conduction and 
the valence band does not influence the spin polarization. The plots show a clear 
polarization around q y = of either +1/2 or —1/2 until the first anti-crossing with a 
higher mode band. 

4. Summary and discussion 

We have studied graphene nanoribbons with armchair boundaries and finite Rashba 
spin-orbit interaction. The solution of the stationary Schrodingcr equation requires the 
choice of a basis set which fulfills the boundary conditions set by the ribbon geometry. 
Therefore the momentum representation of the Hamiltonian is obtained from the 
real space tight binding Hamiltonian via a sine transform, equation (11) and (12), 
instead of the (standard) Fourier transform. We emphasize that this procedure is 
quite different from the case of a graphene armchair nanoribbon without RSOI [23]. 
There the standard approach is to take a basis set consisting of plain waves, which is 
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(a)X=0.01t (b)>.=0.lt 




Figure 7. The spin polarization derived from the tight binding approach is shown 
for (a) A = O.OOlt and (b) A = O.lt and a width of W = 12a In order to compare 
the energy bands with the spin polarization, in the lower panel the spectra for 
the two cases are shown in a color coding agreeing with the color coding of the 
spin polarizations. The lowest mode bands and their polarization are shown in 
red, the second lowest bands and their polarizations are shown in violet and the 
largest shown bands and their polarizations are shown in magenta. 



obtained via a Fourier transform of the bulk problem and to take linear combinations 
of wave function components belonging to the two different K / K 1 points, cf. equation 
(10). But the latter basis cannot be used to diagonalize the RSOI part of the Dirac 
Hamiltonian (5). Due to the derivative operators contained in the expansion up to 
1st order in the momentum of the RSOI Hamiltonian, the use of this standard basis 
is inconsistent and leads to non-Hcrmitian matrices. 

On the other hand, the boundary conditions are naturally satisfied when using 
the sine transform. We have derived the RSOI Hamiltonian in momentum space which 
allows to study the effects of finite RSOI on the energy spectrum in a mode dependent 
manner. Further analytical approximations were presented for the lowest bands and 
the case of small RSOI. We have analyzed the spectra for a variety of ribbon widths 
and RSOI strengths. The main effect is that the spin degeneracy is lifted. For large 
RSOI the spectrum is strongly influenced by the width of the ribbon, while for small 
RSOI this is not the case. The presence of RSOI can lead to the opening of a gap or 
to the appearance of a additional Dirac point. 
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Appendix A. Long wave length expansion of the RSOI Hamiltonian 

The Rashba Hamiltonian in its tight binding representation was presented in (3). We 
use the Fourier transformation of the operators 

a a (Ri) = J d 2 q a a {qy<* R * (A.l) 
bt (Hi + di) = J d 2 q bUq) C -^ R *+ d ^ (A.2) 
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and the relation 



to write 'Hso as: 



s x dj ) • z Sxdy j Sydx.j 



d y ,j ~f~ io] x ,j 



d y ,j ld x ,j 



(A.3) 



Hso =iA / d 2 g^ !,/„,, iad Xtj )e "'' 6^(9)^(9) + H.c. (A.4) 

^ cr 

where er =tl= H — • We use the explicit expressions (1) for the lattice vectors dj and 
after some algebra obtain 



e ^3 



(2 cos(aq x /2 + o-2tt/3> 



+1) + H.C 



(A.5) 



This is the RSOI Hamiltonian as derived in reference [17]. In order to obtain a low 
energy approximation we expand around the K and K' points. The Taylor expansion 
up to first order around the K point {q x = k x — K x and q y = k y ) with a cc = 1 gives 



H 



so 



A / d 2 k< 



(ik x + ky))bUk + K)a t {k + K) + H.c. 



-3iaUk + K)b t {k + K) + R.c. +0 2 (k-K) 



(A.6) 



At the K' point (k x = q x — K' x and k y = q y ) we obtain in a similar fashion: 



n/K' 



\ J d 2 k\ ^(-ik x + k y ))b\(k + K')a±(k + K') + H.c. 



+ 



(-3i)al(fc + K')b±{k + K 1 ) + H.c. + 2 (k - K') 



(A.7) 
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